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We analyze refraction of electromagnetic wave packets on passing from an isotropic positive to 
an isotropic negative refractive index medium. We definitively show that in all cases the energy is 
always refracted negatively. For localized wave packets, the group refraction is also always negative. 



PACS numbers: 78.20.Ci, 41.20.Jb, 42.25.Bs, 84.40.-x 

The existence of a medium with a negative (n < 0) 
index of refraction (NIM), raised several years ago 0, 
has been demonstrated experimentally recently [2]. One 
of the most striking properties of NIM's is that of neg- 
ative refraction for plane waves across the interface be- 
tween positive index materials (PIM) and NIM. Negative 
refraction means that when radiation passes through an 
interface between a PIM and an NIM, the refracted beam 
is on the same side of the normal as the incident beam 
(see Fig. 1), in contrast to the usual positive refraction 
in which they are both the opposit sides of the normal. 

In studies of negative refraction, it is essential to rep- 
resent incident waves as localized wave packets, rather 
than plane waves, since all physical sources of electro- 
magnetic waves produce radiation fields of finite spatial 
and temporal extent because the sources are always of 
finite spatial extent and because they only radiate for 
a finite time. Hence treatments of this problem which 
study waves that extend over infinite distance in all or 
some directions cannot be trusted to reliably predict the 
direction in which a wave will be refracted, and in fact 
treatments based on such extended waves have led to 
a direction of refraction opposite that which one finds for 
spatially localized wave packets, resulting in a great deal 
of controversy and confusion. Although several treat- 
ments using waves of infinite extent in some direction 
(e.g., a plane wave front have obtained negative re- 
fraction, since such a model is unphysical, for the reasons 
given above, we cannot have confidence in conclusions 
obtained from it. 

In this article, we treat refraction of a localized wave 
packet at a PIM-NIM interface both analytically and 
by simulations, demonstrating that it refracts negatively. 
We also present both analytic and numerical studies of 
wave packets constructed from a small number of plane 
waves, on the basis of which we are able to give a plausi- 
ble explanation for why the two plane wave model studied 
by Valanju, et. al.Q|, gives the wrong answer. We find 
that in all cases, including the model of Valanju, et. al., 
the energy and momentum of the wave refract negatively. 
Since electromagnetic waves are detected only when they 
either give up energy to or exert a force on a detector, 
the relevant direction of propagation to consider is that 
of the region of space in which the energy and momentum 
of the wave are nonzero. 



Without sources, Maxwell's equations are V • D = 0, 
V x H = d t T>, V x E + d t B = 0, V ■ B = 0. For plane 
waves of wave vector k and and frequency u>, only three 
equations are independent. Using the usual relationships 
between D(t) and E(t) and between B(t) and H(t)0 
one obtains for such plane waves k x H = — we(o>)E, 
kxE = Ld(j,((j)H. Combining these equations gives us 
a functional relationship between u> and k. Wave propa- 
gation is only permitted for (e, fx, n > 0) or (e, fi, n < 0) 
yjj. In the later case, (E, H, k) will form a left-handed 
triplet while for an ordinary material, (E, H, k) form a 
right-handed triplet. 

A wave packet localized in a compact region of space, 
as occurs in all experimental situations, can be con- 
structed from a continuous distribution of wavevectors. 
Consider such a wave packet incident from outside the 
NIM, E = yE Q j d 2 fc/(k - k )e l ( k - r -"( k )*) with w(k) = 
ck. Here we only consider S'-polarized waves. The 
P-polarized waves can be treated similarly, however. 
Throughout the paper, we choose the z-axis from PIM 
to NIM normal to the interface and the x-axis along the 
interface. If /(k — ko) drops off rapidly as k moves 
away from ko, w(k) can be expanded in a Taylor se- 
ries to first order in k — ko to a good approximation. 
This gives, E = y^oe l(k °' r " w(ko) *g(r - cik /fc ), with 
5 (R) = / d 2 kf{k - k )e i ( k - k °)- R . 

Inside the NIM, k and ko in the argument of the ex- 
ponent get replaced by k r and k r o which are related to 
k and ko by Snell's law 

Kx = k x , k rz = -y/ (n r u/c) 2 - k%. (1) 

Here n r is the refractive index for the NIM and is a func- 
tion of to. Then the wave packet once it enters the NIM 
is given by 

E r = y£oe 4 ^°- r -^ k ^ V(r - v gr t), (2) 

where g r (R) = J d 2 kf(k- k )t k e iR(k ''~ kro) and where 
<k is the transmission amplitude for an incident plane 
wave of wave vector k. It is the standard expression 
for this quantity for the two polarizations of the inci- 
dent plane wavejfj except that the index of refraction 
in each medium is replaced by the index of refraction 
divided by the permeability of the medium. The lat- 
ter modification of the expressions follows quite sim- 
ply from the arguments given in that reference, if one 
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does not replace the ratio of the permeabilities by 1, as 
was done in Ref. 6. Here k r0 denotes k r evaluated at 
k = ko and v gr = Vk r w(k r ) evaluated at k r = k r o. 
Let us expand k r — k r o in the exponential function in 
the expression for g r (R) in a Taylor series in k — ko to 
first order, k r - k r0 s» (k - k ) • Vk(k r - k r0 )|k=k - 
Substituting this in the expression for g r (r), we obtain 
g r (R) = Jd 2 kf(k- k )i k e 4R -[( k - ko > Vk ( k '- k '°)]. If the 
width of the distribution of wave vectors /(k— ko)is small 
compared to the range of k over which ^various signifi- 
cantly, we can to a good approximation simply evaluate 
this quantity at k = ko and put it outside the integral 
over k. Then, the transmission coefficient of the wave 
packet is simply given by |ik | 2 - 

If we carry out the expansion of cj(k r0 )to second order 
in k — k , we are able to show that the wave packet 
spreads out, but if the length and width of the packet 
are much larger than the wavelength corresponding to the 
wave vector ko at the the peak in /(k — ko), we find that 
the amount that the packet spreads out in a given time 
interval is much smaller than the distance traveled by the 
packet in that time. Then clearly under such reasonable 
conditions, the wave packet will remain sufficiently well- 
defined to be able to observe the refraction of the packet. 
The expansion of the frequency in a Taylor series is valid 
for a sufficiently narrow distribution /(k — ko). 

In order to get an explicit expression for g(R), 
let the wave packet have a Gaussian form f(k) = 
(AxAz/tt) exp[— k 2 (Ax) 2 — k 2 (Az) 2 }. Expanding k r in 
a Taylor series around k r o, we get 



g r (R) = cxp[-C 2 x /4(Ax) 2 - C 2 /4(Az) 2 } 



(3) 



with C x = R x + (cn r /v r - l)(k x0 /k rz o)R z , C z = 
(cn r /v r )(k z o/k rz o)R z , and v r = c(dn r uj / 'duj)^ 1 . From 
the above expressions, one can see that the Gaussian 
wave packet moves with v gr . Due to the dispersion, the 
wave packet is deformed in the NIM. 

An NIM is dispersive and causality demands that 
d(ew)/do; > 1 and d(fj,uj) / duj > 1 for nearly transparent 
media 0, Q. For an isotropic NIM, since n r is a func- 
tion of lo only, v gr = c(dn r uj/dui)~ 1 (ck r /n r Lu) = — v r k. r 
with k r the unit vector in the direction of k r . Since v r 
is always positive for transparent media as required by 
causality, the group velocity will be refracted opposite 
the direction of wavevector k r . 

The magnetic field obtained from the electrical field 
through H = (l/w^)k x E is 

Hr = / d 2 kf(k - fc )^4S(k r x y)e * k "— (*>'. 

c J fi r (k) 

from which we find the Poynting vector to be 



E 2 



■ 1 d 2 k / d 2 fc7(k-k )/(k'-k )k r 



x nr ^f\ cos [k r — u}(k)t] cos [k' r ■ r - w(k')t] , (5) 



where we have used the fact that k r ■ y — 0. While there 
is no question that the Poynting vector at a point in a 
medium gives the local direction of energy flow, it does 
not give us the direction of energy flow by a wave packet 
or a group of plane waves as a whole since the direction 
of the Poynting vector varies with space. The integral 
of the Poynting vector over all space, P r = J S r dr, how- 
ever, gives the total momentum carried by a wave packet. 
This quantity divided by the volume over which the wave 
packet is nonzero is the average of the Poynting vector 
over the whole wave packet. Either way, this integral 
clearly represents the direction of motion of the wave 
packet in the medium. From the above expression for 
S r , one has 



(E 2 /2c) / d 2 fc/(k-k ) 2 k 



n r (\k\) 
fi r (\k\)' 



(G) 



Let us consider a coordinate system whose z-axis is along 
ko. The function /(k — ko) 2 will then be a function of k x 
and k z symmetrically peaked around k x = and k z = k^. 
Then writing Eq. (6) as 



P = -§ / ( /- > /,/tk-k„r 



fc^x + k z z 7J r (|fc|) 



we can see that since k is an even function of k x , the 
integrand is an odd function of k x and hence the x- 
component vanishes. Therefore, P, which as argued 
above represents the propagation direction of the wave 
packet, is opposite in direction to ko, i.e., in the direc- 
tion of the group velocity. Hence, the energy refracts 
negatively. 

The negative refraction of the wave packet is illustrated 
by numerical simulation in Fig^ We use the following 
dispersion relation 



n» = -(l/u;)^ 2 - u 2 b )(u* - uft/iu* - u, 2 ). (7) 

for the NIM with ujq < oj < uib- The permeability is 
fx r = (oj 2 — lu 2 )/(uj 2 — u} 2 ). The numbers we used in the 
calculation are, u>o = 1, w& = 3, u p = y/lO, c = 1. FigQ] 
shows stroboscopic snapshots of the electric field inten- 
sity of a propagating wave packet incident on a PIM-NIM 
interface. The negative refraction of the wave packet is 
clearly evident. 

A beam can be constructed as follows: 



E = E J dfc_Le l(ko+ki) - r /(fc±)- 



(8) 



Here kj_ is perpendicular to ko and f(k±) assumes a 
Gaussian form. Note that this construction is different 
from that of Kong et al and Smith et al in that the 
width of the incident packet is made finite in directions 
perpendicular to the direction of propagation. Because 
the NIM is highly dispersive, the incident beam once it 
enters the NIM will no longer be a beam. It will be a 
localized wave packet instead. The electric field E of the 
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PIM NIM 
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FIG. 1: Time-lapse snap shots of the electric field intensity 
of a propagating Gaussian wave packet refracting negatively 
at a PIM-NIM interface. The center wave number is ko — \/E 
with incident angle n/6. The spatial extent of the incident 
wave packet is Aa; = Az — 10. The time step is 50 with speed 
of light c — 1. The dispersion Eq. J7J were used for NIM and 
n = n = 1 for PIM. 
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FIG. 2: Electric field SRi? of a beam with ko = VE and the 
Gaussian weight f(k±) = e _ ' 10fci ' . The incident angle of 
the beam is = 7r/12. 

beam is shown in Fig|2 Just as for the wave packet, the 
beam intensity also refracts negatively. 

We next consider the refraction of wave packets made 
up of a finite number of plane waves. For the cases of 
2 and 3 plane waves analytical expressions are obtained 
for the Poynting vector, momentum and group velocities. 
First consider the case of two plane waves in the xz- 
plane incident from PIM to NIM where the interface is 
at z = 0. Let wave vectors and frequencies be (ki,wi) 
and (ka,W2). We set the polarization in the y-direction. 
Suppose Auj = lo 2 — > 0. The incident wave in PIM 
is E = 2E e lK r - tQt cos(Ak ■ r/2 - Aut/2) with K = 
(ki + k 2 )/2, Q, = (ujx + uj 2 )/2, and Ak = k 2 ki and 
where Eq is the wave amplitude, The electric field of the 
refracted waves is 

E r = 2E' e iKr - r - iQt cos(Ak r • r/2 - Aujt/2)y (9) 



with K r = (k r i +k r 2)/2 and Ak r = k r2 — k r i, where E' Q 
is the amplitude of the electric field and where k r i and 
k r2 are related to ki and k 2 respectively by Eq. (1). 

The relatively long wavelength cosine function in Eq. 
(PJ moves in the NIM with a velocity 

v gr = {Auj/\Ak r \ 2 )(Ak x x + Ak rz z), (10) 

assuming that | Ak| << |K|. From the above expression, 
it is evident that v grx > if Ak x > 0. Since io\ < u>2, 
we have < n(u>x) < n(uj 2 ) and n r (u!i) < n r (uj 2 ) < by 
the requirement of causality which requires d(nuj)/duj > 
0. One has k r -^ z ^r 2 z ~~ ^r 2 ^ 2 x ^ ^* 

Since krz — \krz I, v grz > 0. The group refraction is 
indeed positive [4J. This is due to the simple fact that 
Vgrx > if v gx > 0. Proper dispersion will only give 
v grz > since the energy should propagate away from 
the interface. But we shall see that the above picture is 
not true for the energy flow. 

Let us determine the average Poynting vector (S r ). Us- 
ing the magnetic field corresponding to E r of Eq. (8), 
H r - £oE 3 Li(^^-^z)e !k -- r -^ f / W3 , < S r > is 
found to be given by 

(S r ) - ~(1 + cos A0 r ) |i?g | 2 £ (x^ + i^L) (11) 

3=1 ] 3 

where A<fi r = Ak r ■ r — Aujt. Since k rz < 0, one has 
S x < and S z > 0. Thus, contrary to the refraction 
of the cosine function in Eq. © , the Poynting vector is 
directed in the negative refraction direction, i.e. refracts 
negatively. 

We shall now demonstrate that by including more 
plane waves in our group, one can get negative refrac- 
tion of the group. Actually, just one more plane wave 
can achieve that. Thus, let us include three plane waves, 
whose wave vectors form a triangle, rather than being 
parallel. Let the magnitudes of the wave vectors be k, 
k + Ski, k + Sk2, and their angles with the normal to the 
interface, 9, 9 + S9 1 , 9 + S9 2 . Inside the PIM or the NIM, 
we have 

E = e i(k*x+k z z-wt) fa + gjtp^ _ ct } Skl + ivkSOi) 

+ exp[i(u - ct)5k 2 + ivkS9 2 ]j , (12) 

with u = x sin 9 + z cos 9 and v = x cos 9 — z sin 9 for 
the PIM and u = x sin 9 + az and v = x cos 9 + bz for the 
NIM. Then the lines whose equations are u=constant and 
u=constant are perpendicular for the PIM. Here use has 
been made of the following expansion k rz (k+Sk, 9+59) w 
k rz + o,5k + bk59 with a — /c(sin 2 9 + c|?i r |/u r )/|fc r2 |, b — 
ksm29/2\k rz \- The condition for maximum intensity 
for the quantity in brackets, the long wavelength en- 
velope of the packet, is determined by the equations 
(u — ct)5k\ + vk89\ — 2mi7r, (u — ct)5k 2 + vkS9 2 — 2m 2 TT, 
whose solution in the PIM is x — (ci sin 9 + c 2 cos 9) + 
sin 9 ct,z = (ci cos 9 — c 2 sin 9) + cos 9 ct with c\ = 
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2ir(m 2 89i — 7711662) / (66x5k2 — <5#2<5fci), C2 = 27r(mi<5fc2 — 
m 2 Ski) / (k(S9 1 8k2 — 6628^)), which are clearly only de- 
fined for 5k 1 /5k 2 ^ se 1 /se 2 . 

Inside the NIM, the solution for the location of the 
intensity maxima is x = faa — c\b — bet) /(a cos 9 — 
frsin^), z = (ci cos9~C2 sin 6 + cos 6 ct)/(a cos 9 — b sin 9). 
From the expressions for a and b under Eq. I|12|l. one 
has a, b > and a cos 9 — bs'm9 > 0. Then x(t) and 
z(t), dx/dt < and dz/dt > 0. Thus the refraction will 
be negative. Let the angles of the line u =constant and 
v =constant in the NIM with the z-axis be a and re- 
spectively. Then one has tana = — a/sin#, tan/3 = 
—b/ cos 9 = k x /k rz . So one always has it/2 < a < (3 < ir 
inside the NIM. From the above expressions, one can see 
that the maxima move in the /3 direction, that is, anti- 
parallel to k r . The group velocity in NIM is given by 

V gr = — Vrkr. (13) 

This velocity is independent of how the incident wave 
packet is constructed. The refraction of a group con- 
structed from 4 plane waves is shown in Fig|3| The ar- 
guments presented above demonstrate that for any group 
consisting of 3 or more plane waves whose wavevectors 
are not collinear, the group refraction is negative. 
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FIG. 3: Electric field of negative refraction of 4 plane 
waves with wave vector magnitudes k — Sk,k,k + Sk,k, and 
incident angles, 9,9 — S9, 9,9 + 59, respectively. The center 
wave number is k — v5 with incident angle 9 = 7r/6, 5k = 
0.2, and 59 — n/45. Up to the first order approximation, 
the electric field, Poynting vector, and the moment of this 
group of plane waves are E r = 2e L ^ T (cos vj + cos 5cj> r ), (S r ) = 
— 2(cosci7 + cos 5(j>r) 2 k r /uJ, P£ o11 = — 2Ak. r /cu, respectively. 

While the simulations in Fig[3] clearly show that the 
intensity maxima refract negatively, the normal to the 
planes in which these intensity maxima lie are directed in 



a positive refraction direction. Thus, if one were to imag- 
ine smoothing out all intensity variation in the planes, the 
planes would appear to refract in a positive direction. We 
believe that this is a remnant of the positive refraction 
of the planes of intensity maxima (the cosine function in 
Eq. (PJl) found for the interference patern for the two 
plane wave example of Ref. 3 . When there are only two 
plane waves, this is the only group motion that we see in 
the NIM since for a group consisting of two plane waves, 
there are no intensity variations in these planes. 

One can also look at the energy flow which is repre- 
sented by the Poynting vector. For three wave vectors 
with wave vector magnitudes k — 5k, k, k — 8k, and the 
angles with the normal, 9 — 59, 9 , 9 + 89 respectively, 
the magnetic field can also be calculated and hence the 
Poynting vector up to the first order in both 8k and 89 
is 

(S r ) = — -(1 + 4cos 2 zu + 4costi7COS(50 r )k,./w 

— sinT37sin<5</v (cos#x. + bi)k89/u 

+2 cos 2 w (a - k rz /k)Skz/uj (14) 

where vj = (bz + cos 9 x) k89 and 8(j) r = (az + sin 9 x)8k — 
8wt. Here, (S r ) is not localized; rather it forms a lattice. 
A unit cell is defined as the region in which w changes by 
7T and 8(j) r changes by 27r, as is obvious from the expres- 
sion for E r or (S r ). The area for each unit cell in NIM is 
A = 2iT 2 /(acos9 — bsin9)k8kS9. Instead of integrating 
over all space which will diverge, one can calculate the 
electromagnetic momentum for each cell. Ignoring higher 
order terms in 8k and 89, we get 

P^ c11 = -3A(1 + 2 < 5fc/3fc)K r /2w (15) 

with K r = k r — 2(sin#x + az)8k/3, the average of the 
three wave vectors which make up the group. 

A packet constructed from a finite number of plane 
waves will always give a collection of propagating wave 
pulses with the area of the unit cell inversely proportional 
to 8k and 89. For the above localized waves made of finite 
number of plane waves, the group velocity v gr is parallel 
to P r and anti-parallel to the average wave vector K r . 

In this paper, we have shown that for any localized 
wave packet, the refraction at an interface between a PIM 
and an NIM is always negative. As pointed out earlier, it 
is essential for a correct treatment of this problem to use 
wave packets which are localized in all directions since 
the EM field from any physical source is a localized wave 
packet. 
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